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Abstract. There are several indications that an opaque partonic medium is created in energetic Au+Au 
collisions (^/snn ~ 100 GeV/nucleon) at the Relativistic Heavy Ion Collider (RHIC). At the extreme 
densities of ~ 10 — 100 times normal nuclear density reached even heavy-flavor hadrons are affected sig- 
nificantly. Heavy-quark observables are presented from the parton transport model MPC, focusing on the 
nuclear suppression pattern, azimuthal anisotropy ("elliptic flow"), and azimuthal correlations. Comparison 
with Au + Au data at top RHIC energy ^/sjvjv = 200 GeV indicates significant heavy quark rescatter- 
ing, corresponding roughly five times higher opacities than estimates based on leading-order perturbative 
QCD. We propose measurements of charm-anticharm, e.g., D- meson azimuthal correlations as a sensitive, 
independent probe to corroborate these findings. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 25.75.Ld Collective flow - 25.75.Gz Particle correlations 



1 Introduction 

Recent heavy-ion collisions experiments at the Relativis- 
tic Heavy Ion Collider (RHIC) have generated a lot of ex- 
citement. Among the most remarkable discoveries are the 
large azimuthal momentum anisotropy ( "elliptic flow" ) ^ 
E] and strong attenuation of particles with high transverse 
momentum created in the collision ("jet quenching") [HE], 
which indicate the formation of extremely opaque quark- 
gluon matter that exhibits highly collective, near-hydro- 
dynamic behavior 3!]. The mechanism of rapid random- 
ization and high degree of equilibration in the system 
is not yet understood. For example, it is puzzling that 
dissipative effects (such as viscosity) from nonequilibrim 
transport are significant [Sj and yet ideal (nondissipative) 
hydrodynamics can describe the data quite well[7|. For 
heavy quarks, only partial equilibration is expected be- 
cause collective effects are weaker due to the large mass. 
Therefore, heavy flavor observables are of great interest 
as an orthogonal set of probes to gain more insight and 
cross-check dynamical scenarios. 

There are two main dynamical frameworks to study 
heavy quarks in heavy-ion collisions: parton energy loss 
models |%ll9"lll()llll| and transport approaches [T^lll8lll4llT5l 
Il(3lll7lll8| . Energy loss models consider multiple parton 
scattering in the dense medium in an Eikonal approach 
(i.e., straight-line trajectories), applicable only for very 
large heavy quark energies. The advantage, on the other 
hand, is that coherence effects are taken into account. For 
charm and bottom, small-angle gluon radiation and there- 
fore radiative energy loss is suppressed relative to light 
quarks because of the large quark mass ("dead-cone" ef- 



fect) [HUmi. Surprisingly, recent data from RHIC ggES] 
indicate little light-heavy difference in the high-pT sup- 
pression pattern. Though the puzzle is not resolved yet, it 
is clear now that elastic energy loss, previously neglected, 
plays an important role j 12 1 II 1 1 1 besides radiative energy 
loss. 

Transport models, on the other hand, do not impose 
kinematic limitations but typically include incoherent, elas- 
tic interactions only. They are ideal tools to study equili- 
bration because they have a hydrodynamic (local equilib- 
rium) limit. The dynamics is formulated in terms of (on- 
shell) 6+1D phase space distributions that obey the rel- 
ativitic Boltzmann transport equation |12U18II14II17II18| . 
and the results are mainly sensitive to the transport opac- 
ity of the system ^2]. I n case particles undergo a lot of 
scatterings, the evolution for the bulk of the system (i.e., 
particles that come from "typical" scattering events and 
therefore are affected little by fluctuations in scattering 
angles or the number of scatterings) can be approximated 
with the Fokker-Planck eauation |15IH6| . 

In this work we report on heavy flavor (charm and 
bottom) observables from covariant transport theory with 
elastic 2 — > 2 interactions. The covariant transport solu- 
tions were obtained using the Molnar's Parton Cascade 
(MPC) algorithm 14, 22,. Extending earlier results for charm 
quark elliptic flow (v?) |2H| . we include bottom quarks and 
also study heavy nuclear suppression (Raa) and charm- 
anticharm azimuthal correlations. The results are com- 
pared to RHIC data and also other transport calculations 
in the literature based on the (noncovariant) AMPT trans- 
port model^Z] or the Fokker-Planck limit [16|. 
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2 Covariant transport theory 

We consider here the Lorentz-covariant parton transport 
theory in Refs. 12. 13,24 14,25 , in which the on-shcll par- 
ton phase space densities {fi(x,p)} evolve with elastic 
2^2 rates as 
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+ -Si (a:, Pi) 



(1) 



Here lA^I 2 is the polarization averaged scattering matrix 
element squared, the integrals are shorthands for J 

j d 3 p a /(2E a ), Qi is the number of internal degrees of free- 
dom for species i, while f a> i = fi(x,p a ). The source func- 
tions {Si(x,p)} specify the initial conditions. 

Though, in principle, {*Lj could be generalized for bosons 
and fermions, or inelastic 3 2 processes j26IH8| . no prac- 
tical algorithm yet exists (for opacities at RHIC) to han- 
dle the new nonlinearities such extensions introduce. We 
therefore limit our study to quadratic dependence of the 
collision integral on /. 

We apply JIJ to a system of massless gluons and light 
quarks/antiquarks (q = u, d, s, u, d, s), and charm and bot- 
tom quarks/antiquarks with mass m c = 1.5 GeV, vtih = 
4.75 GeV. All elastic 2^2 QCD processes were taken 
into account: gg — > gg, gQ — > gQ, QQ — > QQ, and 
QQ' — > QQ' . Inelastic 2^2 processes, such as gg «-> QQ, 
are straightforward to include""""") but were ignored here 
for faster simulations. 

The transport solutions were obtained via Molnar 's 
Parton Cascade algorithm |14II22| (MPC), which employs 
the parton subdivision technique |27| to maintain Lorentz 
covariance and causality. Acausal artifacts in the naive 
cascade approach (that uses no subdivision) are known to 
affect basic observables such as spectra, elliptic flow, and 
freezeout distributions in spacetime 13,25 . 

As in Refs. |12U13j . only the most divergent parts of the 
matrix elements were considered, regulated using a Debye 
mass of flu = 0.7 GeV. For perturbative QCD processes at 
leading-order, including scatterings of heavy quarks with 
gluons and light quarks "28], we thus have 
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The last expression was obtained assuming a constant to- 
tal cross section for gg — ► gg (i.e., the weak logarithmic 
energy dependence was neglected). 

In order to reproduce the observed elliptic flow for the 
light parton background at RHIC, scattering cross sections 
between light partons were scaled by a common factor to 
obtain a gg ^gg — 45 mb |13| . about fifteen times the elastic 
2^2 perturbative QCD estimate. This value then fixes 
the total cross sections for all light-parton channels. On 



the other hand, in the spirit of a recent study based on 
Fokker-Planck dvnamics|lfi|. the enhancement of heavy- 
quark cross sections was considered to be a free parame- 
ter. The motivation for this is that these phenomenological 
factors (in part) attempt to incorporate the effect of radia- 
tive processes, which are more important for light partons 
than for the more slowly moving heavy quarks. 

The parton initial conditions for Au + Au at ^snn = 
200 A GeV at RHIC were similar to those in [22], except 
that both initial charm and bottom production was, of 
course, included. For light partons, leading order pQCD 
minijet three-momentum distributions were used (with a 
JT-factor of 2, GRV98LO PDFs, and Q 2 =p 2 T , while Q 2 = s 
for charm). The low-p-r divergence in the jet cross sections 
was regulated via a smooth extrapolation below p± < 2 
GeV to yield a total parton dN(b—0)/dy = 1000 at midra- 
pidity. This choice is motivated by the observed dN c h/dy ~ 
700 and the idea of local parton-hadron duality [2*5). More 
novel hadronization mechanisms, such as parton coales- 
cence, would imply quite different initial conditions |30|. 
Heavy quark momentum distributions were taken from the 
fixed-order plus next-to-leading-log (FONLL) calculation 
in |31| . except for the charm-anticharm correlation results 
in Sec. l3.4l for which correlated c—c distributions were ob- 
tained using the PYTHIA event generator The trans- 
verse density distribution was proportional to the binary 
collision distribution for two Woods-Saxon distributions, 
therefore dNP arton (b=8 ha)/dy » 250. Perfect r? = y cor- 
relation was assumed. 

Because heavy quarks are very rare, scatterings be- 
tween heavy quarks and also the feedback of heavy quarks 
on the light-parton background were neglected. The trans- 
port equations "*TJ then become linear in the heavy quark 
phase space distributions, allowing for weighted test parti- 

cle sampling f(x, p, t) = w l S 3 (x~x l (t))8 3 (p-p l (t)). 

i=l 

The advantage is that sparse, high-p^ phase-space regions 
can be sampled better (the test particle density can be in- 
creased anywhere in phase space, provided the weight is 
reduced in inverse proportion). 



3 Results for heavy flavor 

This section contains heavy flavor results from the trans- 
port model MPC "2*2*| . for conditions expected at RHIC. 
The results below are labelled by the heavy-quark - gluon 
scattering cross section a, for which a wide range was ex- 
plored, from the leading order perturbative QCD estimate 
of a ~ 1.3 mb up to a 15 times enhanced value a = 20 
mb. 



3.1 Nuclear suppression of charm and bottom 

A common observable to characterize parton energy loss 
is the nuclear suppression factor 



Raa(pi 



measured yield in A + A 



expectation for indep. N + N scatterings 



Denes Molnar: Heavy quarks at RHIC from parton transport theory 



3 



which compares the yields to the hypothetical case of in- 
dependent nucleon-nucleon scatterings. In this study, the 
only nuclear affect considered is partonic re-scattering, 
therefore the N + N baseline is given by the initial mo- 
mentum distributions. 

Figure^shows charm and bottom Raa at midrapidity 
as a function of px from covariant transport for Au + Au 
at y/SNN = 200 GeV, with impact parameter 6 = 8 fm. 
At high pt > 5 GeV, heavy-quark yields are suppressed 
because of elastic 2^2 energy loss, and the suppression 
becomes stronger with increasing heavy-quark scattering 
cross section a. At the same a, bottom is less suppressed 
than charm, due to the larger bottom mass. Remarkably, 
already the perturbative a = 1.3 mb generates a signifi- 
cant suppression Raa ~ 0.6 — 0.7. 

At low pt, Raa grows with decreasing px as a nat- 
ural consequence of energy loss. For the largest a = 20 
mb, the charm Raa develops a peak near px ~ 2 GeV, 
which is a clear sign of collective flow (final spectrum has 
a "shoulder-arm" shape due to radial boost). In addition, 
for such large cross sections, the midrapidity (\y\ < 1) 
charm yield is significantly reduced by diffusion in rapid- 
ity, which is the reason why Raa is below one at all px- 

The charm suppression results are qualitatively similar 
to those from the Fokker-Planck approach in JSj. The 
main difference is that charm Raa from the Fokker-Planck 
drops much faster as px increases and does not show any 
sign of charm diffusion in rapidity. This is likely because a 
large final px biases towards fewer scatterings, moreover, 
atypical ("lucky") scatterings contribute significantly to 
the high px yield • These effects reduce the validity of 
the Fokker-Planck approach at high px- 

3.2 Elliptic flow of charm and bottom 

In nonccntral A + A reactions, an independent measure of 
energy loss and deflections in multiple scatterings is dif- 
ferential "elliptic flow", V2{px) = (cos(20)) PT , the second 
Fourier moment of the azimuthal distribution relative to 
the reaction plane at a given px- Figure |21 shows charm 
and bottom V2(px) results at midrapidity for Au + Au at 
^/snn — 200 GeV with 6 — 8 fm from covariant trans- 
port. In the < px < 5 GeV window studied, the re- 
sults are consistent with a monotonic increase with px for 
both charm and bottom. For charm, the increase slows 
down above px > 4 GeV, indicating a turn-over at per- 
haps px ^5 — 8 GeV. At the same px and a, bottom 
i>2 is below charm V2, as generally expected from a mass 
hierarchy observed, in a lower mass region m < 1.5 
GeV, by earlier transport 23 1 7 and ideal hydrodynamic 
calculations jjj . For the perturbative estimate a ~ 1.3 mb, 
charm and bottom elliptic flow are very small, at most a 
few percent. Sizeable heavy-quark elliptic flow V2 ~ 0.1 at 
moderate px ~ 2 — 3 GeV requires 5-10 times enhanced 
cross sections. 

The charm elliptic flow results agree well with ear- 
lier results from covariant transport [23|, and also agree 
within a factor of 2 with results from the AMPT transport 
model |17j . The latter calculation considered quark-quark 



scattering with 3 mb and 10 mb cross sections (no gluons) 
and 2 — 3 times higher parton densities (constituent quarks 
from the "string melting" scenario) , which is roughly equiv- 
alent to the opacities for a ~ 6 — 8 and 20 — 25 mb in our 
case. Charm V2 from AMPT tends to be lower and also 
saturates earlier, around px ~ 2 GeV, whereas our results 
continue to grow until px ~ 3 — 4 GeV. It would be im- 
portant to investigate whether the discrepancy is due to 
differences in initial conditions, or the lack of covariance 
in the AMPT algorithm that has no parton subdivision. 

The results compare qualitatively well to those from 
the Fokker-Planck approach in |16| . The main difference 
is that the Fokker-Planck V2 has a higher slope at low 
Px (i.e., much weaker "mass effect" for charm) and thus 
saturates earlier at high px- In addition, elliptic flow from 
the transport does not exhibit a peak (sharp "rise" and 
"drop") at moderate px, even for the largest cross section 
studied here. 



3.3 Suppression and elliptic flow of decay electrons 

Unfortunately, charm and bottom hadrons are difficult to 
reconstruct experimentally. Though various upgrades are 
under-way to improve detection capabilities, the compro- 
mise at present is to look at "non-photonic" electrons, i.e., 
electrons (and positrons) coming, predominantly, from charm 
and bottom decays. 

Figures and 0] show decay electron results from the 
transport for Au + Au at ^/s^n = 200 GeV. The electron 
(and positron) spectra were calculated via fragmenting the 
heavy quarks into D and B mesons, which were then de- 
cayed using the PYTHIA event generator |32|. Data from 
d + Au collisions at RHIC indicate a very hard heavy- 
quark fragmentation, dominated by momentum fractions 
z ~ 1 35 . For simplicity, we therefore take fragmentation 
functions F c ^d_(z) — 6(1 — z) = Fh^s(z), and consider 
only Z) ± , D°, D° and the corresponding B meson states. 

At high px , electron suppression is very similar in mag- 
nitude to that of heavy quarks, as can be seen in Fig. |31 
The calculations for both 6 = and 8 fm (about 30% cen- 
trality) indicate insufficient suppression for perturbative 
QCD rates. Though experimental uncertainties are large, 
one may speculate that a factor of ~ 5 or higher enhance- 
ment of heavy quark rates is needed for better agreement. 
At low px, much of the structure seen in the heavy quark 
Raa (Fig.[I} gets washed out due to the decay kinematics. 

Figure 0] shows the transport results for the elliptic 
flow of electrons (and positrons) from charm and bot- 
tom decays at midrapidity. Overall, the electron V2(px) 
is very similar to that of charm quarks, except shifted to 
somewhat lower px values (as expected from decays) . This 
corroborates the findings in that only considered elec- 
trons from charm decays. Compared to leading-order per- 
turbative heavy-quark cross sections that give only a few 
percent elliptic flow, a significant V2 ~ 5 — 10% from the 
transport requires a 4 — 8 times increase in heavy-quark 
scattering rates to a ~ 5 — 10 mb. Based on the elec- 
tron Raa data, which suggest a > 5 mb, one expects an 
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Fig. 1. Nuclear suppression factor Raa for charm (left) and bottom quarks (right) as a function of pr in ^4m+^4m at -^/sjvjv = 200 
GeV with 6 = 8 fm, calculated using the covariant transport model MPCJ22] with heavy-quark - gluon scattering cross sections 
a — 1.33 (pluses), 5.33 (open squares), and 20 mb (filled squares). 
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Fig. 2. Differential elliptic flow vi for charm (left) and bottom quarks (right) as a function of pr in Au + Au at y / ijvjv = 200 
GeV with 6 = 8 fm, calculated using the covariant transport model MPG22, with heavy-quark - gluon scattering cross sections 
a — 1.33 (pluses), 5.33 (open squares), and 20 mb (filled squares). 



electron V2 > 5% at moderate pt ~ 1-5 — 5 GeV. Prelim- 
inary data by STAR[37] and PHENIXgHl are compatible 
with an electron elliptic flow of this magnitude but exper- 
imental uncertainties unfortunately prohibit an accurate 
cross-check. 

For the largest a = 20 mb from the transport, at the 
highest pt ~ 4 — 5 GeV the results show a decrease in 
electron elliptic flow. This is because bottom decay con- 
tributions to the overall electron yield start to become sig- 
nificant (eventually take over at higher pr), and bottom 
has a weaker elliptic flow (cf. Fig. [21> ■ 

3.4 Charm-anticharm azimuthal correlations 

Rescatterings not only influence the suppression factor 
and elliptic flow but also the azimuthal correlations be- 
tween two heavy quarks. FigureElshows the charm-antichai 



correlation pattern expected for Au + Au at ^snn — 200 
GeV with b = 8 fm from covariant transport. In the cal- 
culation of this observable, the (correlated) initial charm 
distributions were taken from PYTHIA[3~2"]. PYTHIA pre- 
dicts a strong away-side peak in N + N collisions (i.e., for 
(7 = 0). However, the transport results show a correlation 
strength that is very sensitive to heavy-quark rescatter- 
ing in heavy-ion collisions. The away-side peak is already 
reduced by about half for the small perturbative value 
a ~ 1.3 mb, and as the cross section is increased further, 
the peak rapidly weakens and broadens. Eventually, for 
very large <r, the correlation changes character and a very 
broad near-side peak appears. Measurements of charm- 
anticharm, such as, Z?-meson azimuthal correlations can 
therefore provide an independent way to determine the ef- 
fective heavy quark scattering rates in the dense nuclear 
medium formed in heavy-ion collision. 



Denes Molnar: Heavy quarks at RHIC from parton transport theory 



5 



■J 
ft; 



1.4 
1.2 
1 

0.8 
0.G 
0.4 
0.2 





* PHENIX AuAu. 20-40% 
+ 1.33 mb 
□ 5.33 mb 

■ . 20 mb 

electrons 

+ + + 



Pt [GeV] 



10 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




hh PHENIX AuAu, 0-10% 
+ 1.33 mb 



electrons, 6 = 



10 



Pt [GeV] 



Fig. 3. Left plot shows nuclear suppression factor Raa for decay electrons from charm and bottom as a function of pt in 
Au + Au at y/SNN = 200 GeV with 6 = 8 fm, calculated using the covariant transport model Ml'( ' r 'I'T; with heavy-quark - gluon 
scattering cross sections a — 1.33 (pluses), 5.33 (open squares), and 20 mb (filled squares). The right plot shows the same but 
for b — fm and a — 1.33 mb only. Data on "non-photonic electrons" from PHENIX 19 (filled circles) with statistical errors 
are also shown. 
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Fig. 4. Differential elliptic flow V2 as a function of pr for decay 
electrons from charm and bottom in Au + Au at ^/snn ~ 200 
GeV with b — 8 fm, calculated using the covariant transport 
model MPC|2"2"] with heavy-quark - gluon scattering cross sec- 
tions a — 1.33 (pluses), 5.33 (open squares), and 20 mb (filled 
squares) . 

4 Conclusions 
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Fig. 5. Charm-anticharm azimuthal correlation in Au + Au at 
^/sNN = 200 GeV with 6 = 8, calculated from the covariant 
transport model MPC[2"%] with heavy-quark - gluon scattering 
cross sections a = 1.33 (pluses), 5.33 (open squares), and 20 
mb (filled squares). Charm quarks/anti-quarks with pr > 4 
GeV (triggers) were correlated with (associated) charm anti- 
quarks/quarks with 2 < pr < 4 GeV. The correlation without 
charm rescatterings is also shown (crosses). 



In this work we present heavy-flavor observables in Au + 
Au at RHIC (mainly for impact parameter b = 8 fm) 
from covariant parton transport theory. The heavy-quark 
phase space evolution was studied in an order of mag- 
nitude more opaque light parton (quark and gluon) sys- 
tem than a perturbative parton gas with leading-order 
2-^2 interactions The calculation was driven by the 
cross section a of heavy-quark interactions with gluons. 
The transport solutions were obtained using the covariant 
MPC algorithm[H|E21- 



We find significant charm and bottom suppression with 
Raa ~ 0.5 — 0.65 at high pt > 4 GeV already for leading- 
order heavy-quark matrix elements, and a decreasing Raa 
as the cross section increases (Fig.^l. Electrons from heavy- 
quark decays show a suppression pattern very similar in 
magnitude to that of the heavy quarks (Fig. I3J). Con- 
sistency with data from RHIC 19,20 on "non-photonic" 
electrons requires a roughly five-fold increase in heavy- 
quark opacities relative to perturbative 2 — > 2 scattering 
(more precise data would give better constraints). 
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Charm and bottom elliptic flow V2(pt) from covariant 
transport is at most a few percent for the perturbative es- 
timate of a ~ 1.3 mb. Significant heavy quark elliptic flow 
v 2 > 0.1 at moderate pr ~ 2 — 3 GeV requires about five- 
fold or more enhanced heavy quark rescattering (Fig. |3J) . 
Electron vi (pt) is very similar to that of charm, at least up 
to pt ~ 5 GeV, where bottom contributions to the decay 
electron yield start to become important. For heavy-quark 
opacities indicated by the non-photonic electron Raa data 
at RHIC, the transport predicts v 2 ~ 5 - 10% (Fig. 0J. 
This is within the large uncertainties of current measure- 
ments by STARjHZj and PHENIXPE] - more accurate data 
are highly desirable. 

In addition, we propose a unique observable, charm- 
anticharm azimuthal correlations, as an independent, sen- 
sitive probe of the degree of charm rescatterings in the 
dense parton medium (Fig.^J. Highf>x > 4 GeV charm/anti- 
charm quarks (triggers) were correlated with moderate 
2 < pt < 4 GeV anticharm/charm quarks, both at midra- 
pidity. The strong away-side correlation peak in this ob- 
servable predicted by PYTHIA for N + N collisions is 
strongly reduced (and also broadened) due to rescatter- 
ings, by already a factor of two for the small perturbative 
cross section u ~ 1.3 mb. At very large cross sections, 
the correlation pattern even changes to a broad near-side 
peak. 

We emphasize that this study is limited to 2 — > 2 trans- 
port. Contributions from radiative channels are likely im- 
portant and should be included in the future. The results, 
nevertheless, can serve as a baseline calculation of elastic 
energy loss effects. 

Finally, it would be interesting to extend this calcu- 
lation with hidden heavy- flavor observables, such as J jt\> 
suppression, for which data are also available from RHIC. 
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